16S and ITS data import

Removed taxa not assigned to a phylum. After removing these taxa, 1 16S samples was removed due to having total counts <= 1500.

phyloseq-class experiment-level object
otu_table()   OTU Table:          [ 946 taxa and 48 samples ]:
sample_data() Sample Data:        [ 48 samples by 10 sample variables ]:
tax_table()   Taxonomy Table:     [ 946 taxa by 7 taxonomic ranks ]:
phy_tree()    Phylogenetic Tree:  [ 946 tips and 919 internal nodes ]:
taxa are rows
sample_data() Sample Data:        [ 50 samples by 9 sample variables ]:tax_table()   Taxonomy Table:     [ 246 taxa by 7 taxonomic ranks ]:

Fraction of ASVs classified at each rank

Taxonomic rank	16S	ITS
Kingdom	1.00	1.00
Phylum	1.00	1.00
Class	1.00	0.99
Order	1.00	0.93
Family	0.99	0.89
Genus	0.91	0.86
Species	0.28	0.73

Rarefying to 1569 counts.

$`16S`

otu_table()   OTU Table:          [ 579 taxa and 48 samples ]:
sample_data() Sample Data:        [ 48 samples by 10 sample variables ]:




phyloseq-class experiment-level object
sample_data() Sample Data:        [ 50 samples by 9 sample variables ]:
phy_tree()    Phylogenetic Tree:  [ 111 tips and 109 internal nodes ]:taxa are rows

Predominant phyla

Warning: Returning more (or less) than 1 row per `summarise()` group was deprecated in dplyr 1.1.0.
Please use `reframe()` instead.
When switching from `summarise()` to `reframe()`, remember that `reframe()` always returns an ungrouped data frame and adjust accordingly.Warning: Returning more (or less) than 1 row per `summarise()` group was deprecated in dplyr 1.1.0.
Please use `reframe()` instead.
When switching from `summarise()` to `reframe()`, remember that `reframe()` always returns an ungrouped data frame and adjust accordingly.

Phylum	Total Relative Abundance	Relative abundance B73	Relative abundance CML322
Proteobacteria	0.866	0.829	0.899
Firmicutes	0.065	0.118	0.016
Actinobacteriota	0.038	0.041	0.035
Bacteroidota	0.030	0.010	0.050

Phylum	Total Relative Abundance	Relative abundance B73	Relative abundance CML322
Ascomycota	0.999	0.999	0.999
Basidiomycota	0.001	0.001	0.001

Removed taxa that are present in less than 5% of samples for ASV level dataset. This will be used for differential abundance testing at ASV level.

$`16S`
phyloseq-class experiment-level object
otu_table()   OTU Table:          [ 287 taxa and 48 samples ]:
sample_data() Sample Data:        [ 48 samples by 10 sample variables ]:
tax_table()   Taxonomy Table:     [ 287 taxa by 7 taxonomic ranks ]:phy_tree()    Phylogenetic Tree:  [ 287 tips and 278 internal nodes ]:taxa are rows$ITS

otu_table()   OTU Table:          [ 51 taxa and 50 samples ]:sample_data() Sample Data:        [ 50 samples by 9 sample variables ]:
tax_table()   Taxonomy Table:     [ 51 taxa by 7 taxonomic ranks ]:phy_tree()    Phylogenetic Tree:  [ 51 tips and 50 internal nodes ]:taxa are rows

Agglomerated counts at both genus level and species level.

Warning: Returning more (or less) than 1 row per `summarise()` group was deprecated in dplyr 1.1.0.
Please use `reframe()` instead.
When switching from `summarise()` to `reframe()`, remember that `reframe()` always returns an ungrouped data frame and adjust accordingly.Warning: Returning more (or less) than 1 row per `summarise()` group was deprecated in dplyr 1.1.0.
Please use `reframe()` instead.
When switching from `summarise()` to `reframe()`, remember that `reframe()` always returns an ungrouped data frame and adjust accordingly.Warning: Returning more (or less) than 1 row per `summarise()` group was deprecated in dplyr 1.1.0.
Please use `reframe()` instead.
When switching from `summarise()` to `reframe()`, remember that `reframe()` always returns an ungrouped data frame and adjust accordingly.Warning: Returning more (or less) than 1 row per `summarise()` group was deprecated in dplyr 1.1.0.
Please use `reframe()` instead.
When switching from `summarise()` to `reframe()`, remember that `reframe()` always returns an ungrouped data frame and adjust accordingly.

Most prevalent genera

Genus	Total Relative Abundance	Relative abundance B73	Relative abundance CML322
Pantoea	0.461	0.419	0.500
Klebsiella	0.094	0.024	0.159
Enterobacter	0.051	0.101	0.005
Carnimonas	0.041	0.085	0.000
Burkholderia-Caballeronia-Paraburkholderia	0.030	0.048	0.013
Serratia	0.029	0.033	0.026
Stenotrophomonas	0.029	0.014	0.043
Lactococcus	0.024	0.048	0.003
Sphingobacterium	0.024	0.007	0.039
Achromobacter	0.023	0.008	0.036
Rosenbergiella	0.022	0.047	0.000
Listeria	0.019	0.040	0.000
Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium	0.017	0.009	0.023
Ochrobactrum	0.017	0.013	0.020
Enterococcus	0.010	0.013	0.007

Genus	Total Relative Abundance	Relative abundance B73	Relative abundance CML322
Aspergillus	0.552	0.497	0.608
Sarocladium	0.252	0.226	0.279
Meyerozyma	0.124	0.148	0.101
Talaromyces	0.040	0.076	0.003
Trichoderma	0.013	0.026	0.001
Penicillium	0.006	0.011	0.001
Curvularia	0.005	0.007	0.003
Alternaria	0.004	0.007	0.002
Ramichloridium	0.001	0.001	0.001
Acremonium	0.000	0.001	0.000
Bipolaris	0.000	0.000	0.000
Candida	0.000	0.000	0.000
Ceriporia	0.000	0.000	0.000
Cladosporium	0.000	0.000	0.000
Cyphellophora	0.000	0.000	0.000
Exserohilum	0.000	0.000	0.000
Fomes	0.000	0.000	0.000
Fusarium	0.000	0.000	0.000
Lecanicillium	0.000	0.000	0.000
Moesziomyces	0.000	0.000	0.000
Myrothecium	0.000	0.001	0.000
Paraconiothyrium	0.000	0.000	0.000
Peniophora	0.000	0.000	0.000
Peziza	0.000	0.001	0.000
Phanerochaete	0.000	0.000	0.000
Phlebia	0.000	0.000	0.000
Pyricularia	0.000	0.000	0.000
Scopuloides	0.000	0.000	0.000
Stereum	0.000	0.001	0.000
Tinctoporellus	0.000	0.000	0.000
Trametes	0.000	0.000	0.000
unidentified	0.000	0.000	0.000

Most prevalent species

Species	Total Relative Abundance	Relative abundance B73	Relative abundance CML322
Pantoea_ananatis	0.193	0.207	0.173
Burkholderia_gladioli	0.152	0.195	0.091
Listeria_grayi	0.070	0.119	0.000
Lactococcus_lactis	0.068	0.077	0.056
Sphingobacterium_siyangense	0.050	0.001	0.120
Lactococcus_garvieae	0.043	0.073	0.000
Sphingobacterium_thalpophilum	0.032	0.000	0.079
Sphingomonas_phyllosphaerae	0.030	0.000	0.073
Acinetobacter_baumannii	0.029	0.000	0.070
Pseudomonas_psychrotolerans	0.026	0.008	0.052
Comamonas_sediminis	0.023	0.000	0.056
Corynebacterium_kroppenstedtii	0.021	0.036	0.000
Flavobacterium_anatoliense	0.019	0.000	0.047
Sphingobacterium_multivorum	0.015	0.012	0.019
Staphylococcus_sciuri	0.015	0.026	0.000
Devosia_riboflavina	0.014	0.000	0.035
uncultured_Tistrella	0.013	0.023	0.000
uncultured_bacterium	0.010	0.014	0.003

Species	Total Relative Abundance	Relative abundance B73	Relative abundance CML322
Aspergillus_flavus	0.567	0.525	0.608
Sarocladium_zeae	0.259	0.234	0.284
Meyerozyma_caribbica	0.135	0.168	0.102
Aspergillus_niger	0.011	0.021	0.000
Talaromyces_purpureogenus	0.010	0.019	0.000

Below are barplots of relative taxon abundances for 16S sequencing with samples grouped according to similarity using the neatmap method.

[[1]]

[[2]]

[[3]]

[[4]]

[[5]]

NA

The next two sets are done with ITS counts but made the same way as above.

Warning: stress is (nearly) zero: you may have insufficient data

[[1]]

[[2]]

[[3]]

[[4]]

[[5]]

NA